double mfn(double A, double deltac, double q, double ftNL, char* name, double Bh, double Bf, double mu);

double ratio(double M,double A,double deltac,double q,double ftNL,double z,double B1,double B2, double mu);
double ratio3h(double M,double A, double deltac, double q, double ftNL, double z);
double ratio4h(double M,double A, double deltac, double q, double ftNL, double z, double Bh);
double ratio5h(double M,double A, double deltac, double q, double ftNL, double z, double Bh);
double ratio3f(double M,double A, double deltac, double q, double ftNL, double z);
double ratio4f(double M,double A, double deltac, double q, double ftNL, double z, double Bf);
double ratio5f(double M,double A, double deltac, double q, double ftNL, double z, double Bf);

double Bfeed(double M, double A, double q, double ftNL, double z);

double M2R(double M, double A, double q, double ftNL);
double M3Rh(double M, double A, double q, double ftNL);
double M4Rh(double M, double A, double q, double ftNL, double Bh);
double M5Rh(double M, double A, double q, double ftNL, double Bh);
double M3Rf(double M, double A, double q, double ftNL);
double M4Rf(double M, double A, double q, double ftNL, double Bf);
double M5Rf(double M, double A, double q, double ftNL, double Bf);

double M2RD(double M, double A, double q, double ftNL);
double M3RDh(double M, double A, double q, double ftNL);
double M4RDh(double M, double A, double q, double ftNL, double Bh);
double M5RDh(double M, double A, double q, double ftNL, double Bh);
double M3RDf(double M, double A, double q, double ftNL);
double M4RDf(double M, double A, double q, double ftNL, double Bf);
double M5RDf(double M, double A, double q, double ftNL, double Bf);

double nuc(double deltac, double M, double A, double q, double ftNL, double z);
double nucD(double deltac, double M, double A, double q, double ftNL, double z);

double Bn(double n); double An(double n);

double d2R1(double M, double A, double q, double ftNL);
double d2R2(double M, double A, double q, double ftNL);
double d2R(double M, double A, double q, double ftNL);
double I2R1(double M);
double I2R2(double M);

double d3R1(double M, double A, double q, double ftNL);
double d3R2(double M, double A, double q, double ftNL);
double d3R(double M, double A, double q, double ftNL);
double I3R1(double M);
double I3R2(double M);

double I2R1D(double M);double I2R2D(double M);
double I3R1D(double M);double I3R2D(double M);

double He1(double x); double He2(double x);double He3(double x); 
double He4(double x);double He5(double x); double He6(double x);
double He7(double x); double He8(double x); double He9(double x);
double tgamma(double x);

double const AC = 8.413e-10;
double const h = 0.7;
double const c = 299792.458;
double const Omm = 0.27;
int FEEDER;
int NDFall;

